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1.0  INTRODUCTION 

The  purpose  of  this  project  was  to  determine  the  extent  to  which  quantitative 
information  on  Kelvin  wakes  can  be  obtained  from  synthetic  aperture  radar  (SAR) 
images.  The  method  used  for  extracting  this  information  was  to  model  the  wave 
imaging  process  as  a  linear  system,  as  described  in  the  following  section.  The  wave 
spectrum  is  then  estimated  by  making  corrections  for  the  SAR  ocean  wave  modulatiDn 
transfer  function  and  for  coherent  speckle  effects.  The  methodology  was  applied  to 
both  the  Kelvin  wakes  and  the  ambient  v"  d  waves  observed  during  the  1989  ONR 
Surface  Ship  Wake  experiment.  This  data  set  is  briefly  summarized  and  the  subset  of 
images  analyzed  during  this  project  is  described  in  more  detail  in  section  3.  The 
results  are  presented  in  section  4  and  our  conclusions  are  discussed  in  section  5  of  this 
report. 
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2.0  ANALYSIS  METHODOLOGY 

The  relationship  between  the  SAR  image  intensity  and  the  conditions  at  the  ocean 
surface  is  not  necessarily  unique,  because  of  the  possibility  of  nonlinearitics  in  the 
imaging  process.  However,  under  a  limited  range  of  conditions,  the  imaging  process 
can  be  considered  as  a  linear  or  quasi-linear  system  and  under  these  conditions  the 
system  response  can  be  described  by  means  of  a  modulation  transfer  function  as 
discussed  in  section  2.1  below. 

Variations  in  image  intensity  aie  dso  caused  by  coherent  speckle  effects,  and 
corrections  for  these  effects  must  be  made  in  order  to  extract  quantitative  information 
from  the  images.  These  effects  are  discussed  in  section  2.2. 

Tue  main  focus  of  this  investigation  was  on  the  e..traction  of  information  about  the 
Kelvin  wakes  produced  by  surface  ships.  The  general  characteristics  of  these  wakes, 
and  the  specific  methods  used  to  extract  information  about  them,  are  discussed  in 
section  2.3.  The  methodology  was  applied  to  the  data  set  described  in  section  3,  and 
the  results  are  presented  in  section  4. 

2.1.  LINEARIZED  IMAGING  MODEL 

An  analytical  expression  for  the  Fourier  transform  of  the  SAR  image  intensity  was 
given  by  Lyzenga  (1988)  in  terms  of  the  radar  cross  section  (To(x,y)  and  the  temporal 
correlation  function  p(x,y,r)  of  the  surface  reflectivityl  Assuming  a  nominal 
processor  focus  setting,  this  expression  reduces  to 
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< 


fa/x,y)p(x,y,x)e  '^'^^*'^^'^dxdy 


(1) 


where  r=RK,/2kV,  R  is  the  range  distance,  k  is  the  radar  wavenumber  and  V  i;  the 
platform  velocity.  The  radar  cross  section  and  the  reflectivity  correlation  function  in 
this  case  are  average  values  over  the  SAR  integration  time.  The  reflectivity 
correlation  function  can  be  further  modeled  as 


p(^.y.t)=  e 


(2) 


where  is  the  scene  coherence  time  and  V,(x,y)  is  the  mean  radial  velocity  of  the 
surface  v/ithin  the  resolution  cell  located  at  the  along-track  coordinate  x  and  the 
across-track  coordinate  y.  The  first  term  in  this  equation  describes  the  azimuth  falloff 
effect,  which  limits  the  range  of  along-track  wavenumbers  which  can  be  imaged  due 
to  the  resolution  degradation  caused  by  random  surface  motions.  The  second  term 
leads  to  the  velocity  bunching  effect,  which  enables  waves  to  be  imaged  because  of 
the  orbital  motions  associated  with  these  waves. 

When  either  the  K/V  ratio,  the  azimuth  wavenumber  K„  or  the  radial  velocity  V, 
is  small  enough  so  that  2kV,r  =  RK,Vyv  <  <  1,  the  approximation 
exp(j2kV,r) » 1  -l-j2kV,r  can  be  made.  If,  in  addition,  the  variations  in  radar  cross 
section  are  sufficiently  small  we  can  write  <ro(x,y)=<ro[l+f(>:.y)j  where  f(x,y)<  <  1 
and 
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o,(xO')p(J^.y,T)«  aj[l  *Ax,y)v2kV^(x,y)x]e 


(3) 


neglecting  the  second-order  cross-product  terms.  Substituting  this  into  the  first 
equation  we  obtain  the  linearized  expression 


where  R(K„Ky)  is  the  SAR  modulation  transfer  function  defined  below,  and 


where  S,(K,,Ky)  is  the  wave  height  spectrum,  normalized  so  that  its  integral  over  all 
wavenumbers  is  the  height  variance,  and  6(K)  has  the  properties 


jhWK  ^  \  and  6(0)=I/2n 


(6) 


where  L  is  the  scene  length.  The  SAR  modulation  transfer  function  can  be  written  as 
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D  B 

mJC  rin*4>  ~  — fl/T  cos<J)  cos9  +  i(tn^  sm(^+— flJST  cos<}>  sin<t>  sin0) 


where  m|,*4.5  is  the  hydrodynamic  modulation  transfer  function,  K  is  the  magnitude 
of  the  wavenumber,  tl  is  the  wave  frequency,  ^=tan  '(Ky/KJ  is  the  wave  propagation 
diiection,  B  is  the  incidence  angle,  and 


is  the  tilt  modulation  transfer  function,  where  the  (+)  sign  holds  for  vertical 
polarization  and  the  (-)  sign  holds  for  horizontal  polarization  (Lyzenga,  1986). 

2.2.  COHERENT  SPECKLE  EFFECTS 

We  define  the  image  spectrum  as  Si(K„Ky)=(2T/L)^|F,(K„Ky)|^  where  F,(K„Ky) 
is  the  Fourier  transform  of  the  image  intensity,  i.e. 


where  i(x,y)  is  the  complex  image.  Assuming  that  the  real  and  imaginary  parts  of  the 
complex  image  are  gaussian-distributed  random  variables,  and  applying  the  gaussian 
moment  theorem  (Reed,  1962),  the  expected  value  of  the  image  spectrum  can  be 
written  as 
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=  |<F/ir^;>  |^  +  //  |<i{x.y)/  •(x^x'o'+y^  P«  '^^'^^’^dxdydx'dy'  (JQ) 

where  the  first  term  represents  the  contribution  of  the  resolved  wave  field,  as  modeled 
in  the  previous  section,  and  the  second  term  represents  the  contribution  of  coherent 
speckle  effects  to  the  image  spectrum. 

A  necessary  preliminary  step  in  attempting  to  extract  wave  height  information 
from  the  image  spectrum  is  to  remove  the  speckle  contribution,  as  represented  by  the 
second  term  in  the  above  equation.  The  removal  of  this  speckle  contribution  becomes 
especially  important  at  the  edges  of  the  effective  azimuth  bandpass  region  where  the 
SAR  modulation  transfer  function  becomes  very  small  due  to  the  azimuth  falloff 
effect. 

One  method  of  removing  the  speckle  level  from  the  spectra  is  by  means  of  an 
algorithm  first  described  by  Cordey  and  Macklin  (1989).  This  algorithm  is  discussed 
in  the  following  section.  Another  possible  approach  is  to  use  the  tri-spectrum  or 
fourth-order  moment  sequence  of  the  complex  data.  This  algorithm,  which  is 
discussed  in  section  2.2.2,  has  not  to  our  knowledge  been  previously  applied  to  SAR 
data,  but  is  analogous  to  the  technique  used  in  the  NRL  three-frequency  radar 
(Schuler  et  al,  1985)  to  remove  the  bias  term  from  the  spectrum. 

2.2.1.  Cordey  and  Macklin  Speckle  Removal  Method 

A  method  of  speckle  removal  has  been  described  in  a  recent  paper  by  Cordey  and 
Macklin  (1989).  In  this  method,  which  is  based  on  an  image  model  attributed  to 
Oliver  (1986),  the  speckle  contribution  is  estimated  by  computing  the  complex  image 
autocorrelation  function  <i(x,y)  i*(x-t-x’,y-»-y’)>,  squaring  this  function  and  Fourier 
transforming  it.  The  result  is  basically  equivalent  to  the  speckle  term  obtained  in  the 
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analysis  presented  at  the  beginning  of  this  section.  Thus,  subtracting  this  function 
effectively  removes  the  additive  contribution  of  speckle  from  the  image  spectrum. 
However,  as  shown  by  Oliver  and  discussed  by  Cordey  and  Macklin,  this  function 
may  also  be  used  to  estimate  the  rolloff  of  the  spectrum  due  to  the  stationary-scene 
impulse  response  function  of  the  SAR.  Thus  dividing  the  image  spectrum  by  this 
function  also  corrects  for  this  rolloff  effect  (although  not  for  the  azimuth  falloff  effect 
induced  by  random  surface  motions). 

A  flowchart  of  the  speckle  filtering  process  is  shown  in  Figure  2-1  from  Cordey 
and  Macklin  (1989).  The  estimate  of  the  power  spectrum  (the  right  arm  of  the 
flowchart)  is  formed  in  the  usual  manner  by  forming  the  intensity  image  and  taking 
the  Fourier  transform.  The  estimate  of  the  speckle  filter  (the  left  arm  of  the 
flowchart)  is  formed  by  computing  the  normalized  autocorrelation  function  of  the 
complex  image.  This  autocorrelation  function  is  obtained  by  multiplying  the  Fourier 
transform  of  the  complex  image  by  its  conjugate,  and  then  taking  the  inverse 
transform.  The  autocorrelation  function  is  then  normalized,  multiplied  by  its  complex 
conjugate,  and  Fourier  transformed  again  to  give  an  estimate  of  the  speckle  filter. 

The  filtered  power  spectrum  is  formed  by  dividing  the  power  spectrum  by  the 
speckle  filter.  The  expected  value  of  the  speckle  component  of  the  power  spectrum  is 
now  unity.  This  value  can  now  be  subtracted  from  the  filtered  power  spectrum  to 
yield  a  remainder  which  is  proportional  to  the  wave  height  spectrum,  assuming  the 
imaging  process  is  within  the  linear  regime,  as  discussed  in  section  2. 1  above. 

This  speckle  correction  algorithm  was  applied  to  synthetic  data  to 
quantitatively  assess  whether  the  algorithm  functioned  properly.  To  simulate  a 
speckled  image  of  a  known  signal,  a  256  x  256  complex  array  of  data  was  formed 
using  the  equation 
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Figure  2-1 .  Processing  Scheme  Implemented  for  Speckle  Filtering  Using 
Complex  Imagery  (Extracted  From  Cordey  and  Macklin  [1989]) 
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(11) 


where  a=0.5,  k=T/32,  and  ri(x,y)  is  a  complex  random  variable  whose  real  and 
imaginary  parts  are  zero-mean  Gaussian  random  variables  with  unit  variance. 

Figure  2-2  shows  the  magnitude  of  a  cut  through  the  signal  with  no  noise  present, 
i.e.  with  n  =  l.  Figure  2-3  shows  the  power  spectrum  of  the  signal  with  no  noise. 

Note  that  the  horizontal  axis  of  the  figure  is  in  wavenumber  bins,  and  that  the 
spectrum  has  been  quadrant  swapped  so  that  zero  wavenumber  appears  at  the  center 
(wavenumber  bin  129)  of  the  figure.  Figure  2-4  shows  the  magnitude  of  a  cut 
through  the  signal  with  noise  present,  and  Figure  2-5  shows  the  corresponding  power 
spectrum.  The  signal  components  in  Figure  2-5  are  obvious,  but  note  that  even 
though  the  level  of  the  noise  in  the  spectrum  is  fairly  low,  the  amplitudes  of  the  signal 
components  have  been  modified. 

The  speckle  filter  obtained  using  the  method  described  above  is  shown  in 
Figure  2-6.  Since  the  noise  is  uncorrelated  in  this  example,  the  filter  is  constant  with 
frequency  except  in  the  locations  of  the  signal  components.  Figure  2-7  shows  the 
filtered  power  spectrum.  Note  that  the  mean  of  the  noise  floor  is  approximately  0 
dB,  or  unity,  as  expected.  The  mean  of  the  noise  floor  was  estimated  by  finding  the 
mean  of  thelfirst  100  frequeny  bins,  and  this  value  was  subtracted  from  the  filtered 

I 

power  spectrum.  Following  the  procedure  described  by  Rotheram  and  Macklin 
(1985),  the  filtered  power  spectrum  was  then  thresholded  at  a  value  three  times  the 
mean  value,  ii^^e.  if  the  value  of  the  power  spectrum  was  less  than  threshold  value,  the 
value  was  set  i,o  zero.  Figure  2-8  shows  the  thresholded  power  spectrum. 

Finally,  the  thresholded  power  spectrum  was  multiplied  by  the  speckle  filter  to 
recover  the  original  values  of  the  signal  component,  which  are  shown  in  Figure  2-9. 
Comparing  Figure  2-9  with  Figure  2-3,  it  can  be  seen  that  a  power  spectrum 
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Figure  2-4.  Magnitude  of  a  Cut  Through  the  Synthetic  Signal  With  Noise 
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essentially  equivalent  to  the  power  spectrum  of  the  signal  without  noise  is  obtained. 
The  results  of  this  example  suggest  that  the  speckle  correction  method  can  be  used  to 
quantitatively  estimate  the  signal  component  of  the  image  spectrum. 


2.2.2.  Trispectrum  Algorithm 

Another  approach  to  the  problem  of  removing  speckle  effects  was  investigated 
during  this  study,  using  an  alternative  algorithm  for  estimating  the  image  spectrum. 
The  method  bears  some  similarity  to  the  NRL  three-frequency  technique  (Schuler  et 
al,  1985)  in  that  the  expected  value  of  the  spectral  estimate  is  simply  proportional  to 
the  wave  spectrum,  and  has  no  background  or  speckle  contribution.  However,  the 
method  has  some  other  disadvantages  which  appear  to  limit  its  usefulness,  as 
discussed  below. 

Using  the  Fourier  transform  product  or  convolution  theorem,  the  Fourier 
transform  of  the  image  intensity  can  be  written  as 


where  f(K„Kj,)  is  the  Fourier  transform  of  the  complex  image  i(x,y).  Thus,  the 
image  spectrum  can  be  written  as 
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5(/:^P=  /// jAKl^lTiK-KX-K;if\K'Jjc!;)j{K'J-K^'J-K)dK^^^^  .  (13) 

Applying  the  gaussian  moment  theorem  to  the  integrand,  the  exnected  value  of  this 
spectrum  can  be  seen  to  consist  of  a  "signal"  component  and  a  bias  term  as  discussed 
at  the  beginning  of  section  2.2.  However,  if  we  define  the  quantity 


/ lf\Kl-KX-f^y(K'XTi.K^KX*K)  dK'^Kl  (14) 


and  apply  a  similar  analysis,  we  find  that  the  gaussian  moment  decomposition  rf  suits 
in  two  identical  terms,  both  of  which  are  proportional  to  the  "signal"  spectrum  under 
the  assumption  that 


(15) 


which  follows  from  the  above  definitions.  Thus,  the  quantity  SjCK^Ky),  which  we 
refer  to  as  the  "tri-spectrum"  below,  can  be  used  as  an  unbiased  estimator  of  the  wave 
spectrum,  i.e.  no  corrections  for  speckle  effects  need  to  be  applied. 

This  method  was  tested  using  one-dimensional  simulated  data  generated  using  the 
same  model  as  discussed  in  section  2.2.1  but  with  128  points  and  values  of  a=0.2 
and  0.5.  Figures  2-10  and  2-11  show  comparisons  between  the  tri-spectrum  and  the 
conventional  spectrum  for  these  cases,  averaged  over  128  realizations.  The  peak  due 
to  the  sinusoid  appears  in  the  same  spectral  location  for  both  methods  and  also  has 
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about  the  same  amplitude.  The  DC  component  for  both  spectral  forms  has  been 
removed. 

As  these  figures  illustrate,  the  average  background  level  is  non-zero  for  the 
conventional  spectrum  and  essentially  zero  for  the  tri-spectrum.  However,  it  is  also 
clear  from  these  figures  that  the  variance  of  the  background  is  significantly  increased 
for  the  tri-spectrum.  Therefore,  it  does  not  appear  that  the  tri-spectrum  method  has 
any  practical  advantage  over  the  conventional  spectrum  with  an  estimate  of  the  bias 
removed  as  discussed  in  section  2.2.1. 

2.3.  KELVIN  WAKE  MODEL 

The  basic  condition  which  defines  the  stationary  pattern  of  waves  generated  by  a 
moving  ship  is  that 


Vp  =  V,  cos  <t>  (16) 

where  Vp  is  the  phase  velocity  of  any  given  wave  component,  </>  is  the  propagation 
direction  of  this  wave  component  relative  to  the  ship  track,  and  V,  is  the  ship  speed. 
This  relationship  can  also  be  written  in  the  form 

V.-K  =  0  =  fro  (17) 

where  V,  is  the  ship  velocity  vector,  K  is  the  wavenumber  vector,  V,*K  denotes  the 
dot  product  of  these  two  vectors,  and  0  =  l(gK)  is  the  wave  frequency.  The  set  of 
waves  making  up  the  Kelvin  wake,  as  defined  by  either  of  these  equations,  is  shown 
graphically  in  Figure  2-12  for  a  range  of  ship  speeds. 


Range  Wavenumber  (rads/m) 

0.4a  -0.36  -0.24  -0.12  0.00  0.12  0.24  0.36 


Kelvin  Wake  Speclruiii 


Figure  2-1 2.  Schematic  of  Kelvin  Wake  Energy  Locations  in  Spectral  Space 


A  further  analysis  involving  the  group  velocity  of  the  waves  can  be  carried  out 
(e.g.  Tricker,  1964)  in  order  to  determine  the  position  of  each  wave  component  in  the 
wake  pattern.  Such  an  analysis  shows  that  the  waves  which  propagate  furthest  from 
the  track,  in  terms  of  angular  position  relative  to  the  ship,  are  those  having  a 
propagation  direction  of  0= ±35.26°  relative  to  the  ship’s  heading.  Waves  with 
smaller  values  of  <i>  are  referred  to  as  transverse  waves,  and  those  with  larger  values 
of  4>  are  called  the  diverging  components.  At  any  given  instant  of  time,  the  entire 
wake  pattern  is  spatially  confined  to  a  wedge  having  a  half-angle  of  19.47  degrees 
with  its  apex  at  the  ship’s  position. 

The  kinematic  analysis  summarized  above  determines  the  locus  of  points  in 
wavenumber  space  which  make  up  the  Kelvin  wake  for  a  given  ship  speed,  as  shown 
in  Figure  2-12.  However,  it  does  not  specify  the  amplitude  of  these  components.  In 
general,  the  amplitude  distribution  along  the  lines  shown  in  Figure  2-12  depends  on 
the  ship  hull  size  and  shape,  as  well  as  the  ship  speed.  Modeling  of  this  amplitude 
distribution  is  clearly  beyond  the  scope  of  this  report. 
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3.0  DATA  SET  DESCRIPTION 

The  data  analyzed  in  this  report  was  collected  during  a  two-week  ship  wake 
experiment  conducted  off  the  California  coast  in  January  and  February  1989  (see 
Chapman,  1989)  for  a  description  of  this  experiment).  The  experiment  utilized  6 
different  ship  types  and  the  data  was  collected  over  a  wide  range  of  environmental 
conditions.  A  total  of  eight  data  collection  flights  consisting  of  113  passes  were 
conducted  with  the  ERIM/NADC  P-3/SAR  system.  The  P-3/SAR  system  (insists  of 
a  three  frequency  polarimetric  SAR  installed  in  a  Naval  Air  Development  Center 
(NADC)  P-3  aircraft.  The  system  is  a  fully  focused  side-looking  SAR  capable  of 
looking  out  either  side  of  the  aircraft.  The  system  can  operate  simultaneously  at  X- 
,C-,  and  L-band  (9.4,  5.3,  and  1.2  Ghz,  respectively).  The  system  transmits  either 
horizontally  -  (H)  or  vertically-polarized  (V)  pulses  and  can  receive  either  or  both 
polarizations.  The  P-3  was  navigated  during  SAR  data  collection  by  an  Inertial 
Navigation  System  (INS)  updated  by  a  LORAN-C.  A  more  complete  description  of 
the  P-3/SAR  system  is  given  by  Sullivan,  et.  al.  (1987). 

During  the  ship  wake  experiment,  the  P-3  SAR  was  operated  in  the  configuration 
shown  in  Table  1.  Detailed  analyses  were  carried  out  during  the  present  study  using 
four  data  sets  collected  on  28  January  1989,  designated  as  runs  4-1,  6-4,  7-2  and  7-4. 
The  pertinent  aircraft  and  ship  operating  parameters  for  these  runs  are  given  in 
Table  2.  The  ship  used  during  these  runs  was  the  USS  Chandler  (DDG). 

During  run  4-1,  which  was  collected  at  approximately  11:11  PST,  the  winds 
reported  by  the  Chandler  were  about  5  knots  in  the  area  just  south  of  Santa  Cruz 
Island  where  the  data  was  collected.  The  other  runs  were  collected  further  offshore, 
as  shown  in  Figure  3-1,  at  approximately  14:20  PST  (run  6-4),  15:27  PST  (run  7-2) 
and  16:02  PST  (run  7-4).  The  winds  during  this  time  interval  were  fairly  constant,  as 
measured  by  the  R/V  McGaw,  at  about  10  knots  from  300°N  (Gasparovic  and 
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TABLE  1 

P-3/SAR  CONFIGURATION  DURING  SHIP  WAKE  FLIGHTS 


Mode 


Altitude 


Ground  Speed 


Resolution 


Swath  Width 


Range  Samples 


!  Slant-Range  Sampling  Rate 


Along-Track  Sampling  Rate 


Xvv«  Lvv»  Cvv»  I-^ 


10  Kft 


250  Kts  (nominal) 


2.5  m  (range)  x  2.1  m  (azimuth) 


9.8  km  (slant-range) 


4096 


2.4  m 


0.54  m 


TABLE  2 

SAR  AIRCRAFT  AND  SHIP  PARAMETERS 
FOR  DATA  SETS  ANALYZED  IN  THIS  REPORT 
(28  JANUARY  1989) 


RUN/PASS 


A/C 

HEADING 

SHIP 

HEADING 

A/C  SPEED 

SHIP  SPEED 

090 

270 

146.7  m/s 

25  kts 

210 

210 

143.6  m/s 

18  kts 

255 

345 

155.5  m/s 

25  kts 

165 

345 

123.2  m/s 

25  kts 
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Figure  3-1.  Target 
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Drummond,  1989).  Ambient  wave  spectra  were  also  measured  with  a  spar  buoy 
launched  by  the  R/V  McGaw  (Gasparovic,  1989).  The  spectrum  measured  during  the 
time  interval  from  15:05-15:25  is  shown  in  figure  3-2.  This  spectrum,  which  is  very 
similar  to  the  one  measured  during  the  previous  period  from  14:30-14:50,  shows  two 
peaks  of  roughly  the  same  amplitude  at  frequencies  of  0.06  Hz  and  0.125  Hz, 
corresponding  to  wavelengths  of  430  meters  and  100  meters,  respectively.  The 
significant  wave  height  was  approximately  1  meter. 

The  X-band  (VV)  and  L-band  (VV)  images  for  these  three  passes  are  shown  in 
Figures  3-3  through  3-6.  These  images  have  been  artificially  flattened  to  enhance 
feature  visibility.  The  spaital  coverage  of  these  images  is  6  km  in  azimuth  (left  to 
right)  and  10  km  in  slant  range.  Further  processing  and  analysis  of  these  images  is 
discussed  in  the  following  section. 
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Figure  3-2.  Ambient  Wave  Height  Spectrum  Measured  by  APL  Near  the 
R/V  Me  Gaw  (From  Gasparovic,  1989) 
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Pt.  Mugu  January  28,19<39 
DDG  Run  6-4  X-VV 
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Pt.  Mugu  January  28,1989 
DDG  Run  6-4  L-W 


Figure  3-4.  X-VV  and  L-VV  Images  for  Pass  6-4 
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Pt.  Mugu  January  28, 1989 
DDG  Run  7-2  X-VV 
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4.0  DATA  ANALYSIS  AND  RESULTS 

Image  intensity  spectra  were  generated  from  the  SAR  data  sets  described  in  the 
previous  section,  over  areas  containing  the  Kelvin  wake  from  the  USS  Chandler  and 
in  background  areas  containing  only  the  ambient  waves.  Several  types  of  analysis 
were  then  carried  out  to  extract  wave  height  information  from  both  the  wake  and  the 
ambient  waves.  Section  4. 1  describes  an  analysis  performed  on  the  ambient  wave 
spectra  in  order  to  determine  the  adequacy  of  the  linear  imaging  model  discussed  in 
section  2.  This  section  also  describes  the  application  of  the  Cordey-Macklin  speckle 
correction  algorithm  to  this  data  set.  Section  4.2  then  addresses  the  inverse  problem 
of  estimating  wave  amplitudes  from  SAR  imagery  by  correcting  for  the  modulation 
transfer  function. 

4.1.  MODEL  VALIDATION 

The  ambient  wave  data  set  was  used  to  test  the  Cordey-Macklin  speckle 
correction  algorithm,  as  discussed  in  the  following  section.  The  image  spectra  for 
passes  6-4,  7-2  and  7-4  were  then  compared  with  predictions  from  the  linear  imaging 
model  discussed  in  section  2  in  order  to  test  this  model.  This  comparison  is  presented 
in  section  4.1.2  below. 

4.1.1.  Speckle  Correction 

The  speckle  reduction  algorithm  described  in  section  2.2.1.  was  applied  to  a  512 
X  512  pixel  subset  of  the  L-band  VV  polarization  SAR  image  for  pass  6-4  in  order  to 
test  the  algorithm.  The  pixel  spacings  in  the  image  were  2.4  m  in  range  and  1.62  m 
in  azimuth.  The  complex  image  was  first  normalized  such  that  the  mean  of  the 
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intensity  image  was  unity.  The  power  spectaim  of  the  image  was  then  formed  by 
taking  the  squared  magnitude  of  the  complex  image,  Fourier  transforming  this  image 
intensity,  and  taking  the  squared  magnitude  of  tb;  Fourier  transform.  Figure  4-1 
shows  a  cut  through  pnwpr  snectnim  alone  the  azimuth  wavenumber  axis,  and 
Figure  4-2  shows  a  cut  along  the  range  wavenumber  axis.  The  horizontal  axis  in  the 
figures  is  in  wavenumber  bins,  and  the  spectral  quantities  have  been  quadrant 
swapped  suL.h  that  zero  wavenumber  is  in  the  center  of  the  range  (wavenumber  bin 
257).  Note  that  in  both  figures  the  noise  floor  falls  off  with  increasing  wavenumber 
due  to  the  correlation  between  pixels  in  the  image. 

A  two-dimensional  color  representation  of  the  power  spectrum  is  shown  in  the 
upper  left  hand  comer  of  Figure  4-3.  The  range  direction  is  vertical,  and  the  azimuth 
direction  is  horizontal.  The  color  scale  used  in  the  image  (and  in  all  images  that  will 
be  presented  in  this  section)  is  a  60  dB  scale  fixed  at  the  maximum  value  in  the 
image.  The  color  scale  runs  from  black  and  purple  at  the  low  end  of  the  scale  to  red 
at  the  high  end.  Green  represents  the  midpoint  of  the  scale.  The  power  spectrum  in 
Figure  4-3  shows  strong  signal  returns  aligned  in  the  range  direction.  The 
surrounding  noise-dominanted  part  of  the  spectrum  falls  off  radially  with  frequency. 
The  variation  of  the  noise  component  of  the  spectrum  with  frequency  makes  it  hard  to 
estimate  and  remove,  which  is  the  motivation  for  the  algorithm  discussed  here. 

A  speckle  filter  for  the  image  was  formed  in  the  manner  described  in  section 
b.2.1.  Figure  4-4  shows  a  cut  through  the  speckle  filter  along  the  azimuth 
jvavenumber  axis,  and  Figure  4-5  shows  a  cut  along  the  range  wavenumber  axis.  A 
two-dimensional  color  representation  of  the  speck’'*,  filter  is  shown  in  the  upper  right 
h^d  comer  of  Figure  4-3.  Note  that  the  slow  variation  of  the  speckle  filter  with 
frMuency  has  roughly  the  same  shape  as  the  variation  of  the  noise  floor  in  the  power 
spajtnim. 
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Figure  4-2.  A  Cut  Through  the  Collected  SAR  Image  Power  Spectrum  Along  the  Range 
Wavenumber  dc 
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Figure  4-3.  Collected  SAR  l.nage  (Upper  Left)  Power  Spectrum.  (Upper  Right)  Speckle 
Filter,  (Lower  Left)  Filtered  Power  Spectrum,  and  (Lower  Right)  Filtered 
Power  Spectrum  After  Thresholding.  Al!  Data  are  Displayed  on  a  60  dB  Scale 
Fixed  at  tne  Maximum  in  the  Individual  Data  Set.  The  Color  Scale  Runs  From 
Black  and  Purple  at  the  Low  End  to  Red  at  the  High  End.  Green  is  in  tne 
Micdle  of  the  Scale. 
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A  filtered  power  spectrum  was  formed  by  dividing  the  power  spectrum  by  the 
speckle  filter.  Figure  4-6  shows  a  cut  through  the  filtered  power  spectrum  along  the 
azimuth  wavenumber  axis,  and  Figure  4-7  shows  a  cut  along  the  range  wavenumber 
axis.  A  two-dimensional  color  representation  of  the  filtered  spectrum  is  shown  in  the 
lower  left  hand  comer  of  Figure  4-3.  Note  in  these  figures  that  the  noise  floor  is  now 
relatively  constant  with  frequency,  with  the  mean  of  the  noise  floor  approximately 
0  db  (unity).  The  cut  along  the  azimuth  wavenumber  axis  (Figure  4-6)  does  show 
more  variation  with  frequency  than  the  cut  along  the  range  wavenumber  axis 
(Figure  4-7),  suggesting  that  the  filter  works  better  in  range  than  in  azimuth. 

Now  that  the  variation  of  the  noise  floor  with  frequency  is  relatively  constant,  it 
can  be  more  reliably  estimated.  Two  estimates  of  the  noise  floor  were  obtained  by 
finding  the  mean  of  power  spectrum  values  for  azimuth  freouency  bins  less  than  50 
and  less  than  100.  As  discussed  above,  the  expected  value  should  theoretically  be 
unity.  The  average  of  bins  less  than  50  yielded  1.05,  and  the  average  of  bins  less  than 
1(X)  yielded  1.08.  1.08  was  chosen  as  the  mean  and  was  subtracted  from  the  filtered 
power  spectrum. 

The  next  step  in  the  process  is  to  threshold  the  resulting  filter  power  spectrum. 
The  correct  threshold  value  is  a  subject  of  research.  To  investigate  tl:e  effects  of 
changing  the  threshold,  the  filter  power  spectrum  data  were  thresholded  at  values  one, 
three,  five  and  seven  times  the  estimated  mean  value  of  the  noise  floor.  Figure  4-8 
shows  the  resulting  spectra.  The  upper  left  spectrum  was  thresholded  at  the  estimated 
mean  value  of  the  noise  floor,  the  upper  right  spectrum  was  thresholded  at  three  times 
the  estimated  mean  value  of  the  noise  floor,  the  lower  left  spectrum  was  thresholded 
at  five  times  the  estimated  mean  value  of  the  noise  floor,  and  the  lower  right 
spectrum  was  thresholded  at  seven  times  the  estimated  mean  value  of  the  noise  floor. 
Qualitatively,  it  appears  that  thresholding  the  filter  power  spectrum  at  one  times  the 
mean  is  too  low  since  much  of  the  noise  remains.  On  the  other  hand,  thresholding  the 
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Figure  4-6.  A  Cut  Through  the  Collected  SAR  Image  Filtered  Power  Spectrum  Along  the 
Azimuth  Wavenumber 


Figure  4-8.  Filtered  SAR  Irr.age  Power  Spectrum  After  Thresholding  at  (Upper  Left)  One, 
(Upper  Right;  Three,  (Lower  Left)  Five,  and  (Lower  Right)  Seven  Times  the 
Estimated  Mean  Value  of  the  Noise  Floor.  All  Data  are  Displayed  on  a  60  dB 
Scale  Fixed  at  the  Maximum  in  the  Individual  Data  Set.  The  Color  Scale  Runs 
From  Black  and  Purple  at  the  Low  End  to  Red  at  the  High  End.  Green  is  in  the 
Middle  of  the  Scale. 
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filtered  power  spectrum  at  seven  times  the  mean  may  remove  some  signal  component. 
More  work  needs  to  be  done  to  determine  the  optimum  threshold  level. 

4.1.2.  MTF  Analysis 

The  two-dimensional  image  spectra  for  the  ambient  waves  in  passes  6-4,  7-z,  and 
7-4  on  28  January  1989  are  shown  in  Figures  4-9  through  4-11.  These  spectra  were 
obtained  from  image  subsets  centered  at  an  incidence  angle  of  approximately  55 
degrees,  and  a  range  distance  of  about  5.4  km  in  each  image.  The  primary  difference 
from  pass  to  pass  is  the  orientation  of  the  flight  path  with  respect  to  the  direction  of 
the  dominant  wave  pattern,  which  appears  to  be  roughly  coincident  with  the  wind 
direction  of  300“.  Thus,  the  dominant  wave  direction  was  approximately  in  the  range 
direction  during  run  6-4,  and  was  roughly  45  degrees  on  either  side  of  the  range 
direction  during  runs  7-2  and  7-4.  Both  HH  and  VV  imagery  were  collected  during 
all  ibree  runs.  Since  the  ambient  wave  spectrum  was  roughly  the  same  for  all  three 
passes,  the  differences  among  the  image  intensity  spectra  can  be  attributed  to 
variations  in  the  SAR  response  function  with  collection  geometry  and  polarization. 

The  behavior  of  the  image  intensity  spectrum  was  investigated  as  a  function  of 
azimuth  angle  and  polarization.  To  examine  this  behavior,  the  image  intensity 
spectrum  was  integrated  over  wavenumber  and  plotted  as  a  function  of  azimuth  angle. 
The  spectra  were  integrated  in  wavenumber  between  2^/150  rad/m  and  2x/50  rad/m 
in  order  to  encompass  the  peak  wavenumber.  The  integrated  spectra  for  HH  and  VV 
polarization  are  given  in  Figures  4-12  and  4-13  respectively.  As  can  be  .seen  from  the 
figures,  the  integrated  spectra  have  distinctly  different  shapes  foi  different  look 
directions  and  polarizations. 

One  particularly  interesting  feature  of  the  integrated  spectra  is  the  deep  null 
observed  near  the  range  wavenumber  axis  in  run  6-4.  This  effect  was  first  observed 
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Figure  4-9.  L-VV  and  L-HH  Image  Subsets  and  Image  Spectra  Extracted  From  Run  6-4 


Figure  4-10.  L-VV  and  L-HH  Image  Subsets  and  Image  Spectra  Extracted  From  Run  7-2 
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Figure  4-1 1 .  L-VV  and  L-HH  Image  Subsets  and  Image  Spectra  Extracted  From  Run  7-4 
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Figure  4-12.  Empirical  HH  Image  Spectra  From  Runs  6-4,  7-2  and  7-4  Integrated  Over 
Wavenumber 
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Figure  4-13.  Empirical  VV  Image  Spectra  From  Runs  6-4, 7-2  and  7-4  Integrated  Over 
Wavenumber 
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by  Bruning  et  al  (1988)  and  was  attributed  to  an  interference  effect  in  the  SAR 
modulation  transfer  function  (MTF).  Note  that  the  null  in  the  VV  spectrum  is  deeper 
than  that  in  the  HH  spectrum.  Another  observation  is  that  the  peaks  of  the  integrated 
spectra  for  runs  7-2  and  7-4  are  less  than  90°  apart,  indicating  that  these  spectra  are 
rotated  toward  the  range  direction. 

To  investigate  the  angular  variations  of  the  image  intensity  spectra  shown  in 
Figures  4-12  and  4-13,  the  image  intensity  spectrum  was  simulated  by  assuming  a 
simple  form  for  the  wave  height  spectrum  and  multiplying  by  the  SAR  MTF 
discussed  in  section  2.1.  The  same  height  spectrum  was  used  to  simulate  the  intensity 
spectrum  for  all  three  passes,  thus  variations  in  the  shape  of  the  intensity  spectrum 
that  we  calculate  using  this  procedure  are  due  to  variations  in  the  MTF.  From  the 
discussion  in  section  2.1,  we  would  expect  this  procedure  to  be  valid  if 
RK,VyV  <  1/2. 

The  R/V  ratios  for  the  image  subsets  used  to  calculate  the  spectra  in  Figures  4-12 
and  4-13  are  37.8  sec  for  run  6-4,  35.3  sec  for  run  7-2,  and  43.4  sec  for  run  7-4,  the 
variation  being  due  primarily  to  differences  in  the  platform  velocity.  An  estimate  of 
the  maximum  radial  velocity  can  be  obtained  by  noting  that  for  a  Pierson-Moskowitz 
spectrum  the  ratio  of  the  velocity  variance  to  the  height  variance  is 
approximately  60  f,,^  where  f^  is  the  peak  frequency.  Using  fo=0.125  Hz  and 
(r|,=0.25  m,  we  obtain  (r„=0.24  m/sec.  Equating  with  V„  the  parameter 
a=RK,Vr/V  has  the  value  0.46  at  the  peak  of  the  spectrum  for  run  7-4,  which  is  the 
"worst  case"  (largest  R/V).  We  can  therefore  expect  the  linear  imaging  model  to  be 
valid  for  wavenumbers  near  the  spectral  peak,  although  it  is  approaching  the  limits  of 
its  validity  in  some  cases. 

In  order  to  test  the  predictions  of  this  theory,  the  wave  height  spectrum  was 
modeled  as  a  simple  gaussian  function 
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having  a  width  ctk  and  centered  at  the  wavenumbers  K„=K^cos<^o  and  Kyo=KoSin<^o, 
where  Ko  is  the  peak  wavenumber  and  0^  is  the  peak  direction.  The  peak  wavelength 
was  chosen  as  80  meters,  based  on  the  observed  location  of  the  peak  in  the  image 
spectrum,  and  the  other  parameters  were  iterated  to  yield  the  best  fit  to  the  observed 
angular  spectra.  This  iteration  resulted  in  </>o  being  chosen  as  5“  from  the  nominal 
wind  direction,  and  the  spectral  width  being  chosen  as  <^=0.0313  rad/m.  The 
amplitude  A  was  taken  as  1.0  x  10'^  m^.  The  modeled  wave  height  spectra  integrated 
over  wavenumber  are  plotted  in  Figure  4-14  as  a  function  of  the  azimuth  angle 
relative  to  the  flight  direction  for  each  pass. 

To  model  the  image  intensity  spectrum,  the  wave  height  spectrum  was  multiplied 
by  the  squared  magnitude  of  the  SAR  MTF  defined  in  section  2.1.  The  coherence 
time  Te  was  taken  to  be  0.05  sec.  The  resulting  integrated  image  intensity  spectra  for 
HH  and  VV  polarization  are  shown  in  Figures  4-15  and  4-16,  respectively.  Note  that 
qualitatively,  the  behavior  of  the  spectra  as  a  function  of  azimuth  angle  and 
polarization  in  Figures  4-15  and  4-16  is  essentially  the  same  as  the  behavior  of  the 
spectra  in  Figures  4-12  and  4-13.  The  shape  and  relative  amplitude  of  the  modeled 
spectra  show  the  same  variation  with  polarization  and  azimuth  angle  as  the  spectra 
extracted  from  the  actual  imagery.  Most  notably,  the  modeled  spectra  for  Run  6-4 
show  the  null  about  the  range  direction,  with  the  null  in  the  VV  spectrum  being 
deeper  than  the  null  in  the  HH  spectrum. 

The  behavior  of  the  spectra  as  a  function  of  azimuth  angle  and  polarization  can  be 
simply  explained  by  examining  the  HH  and  W  modulation  transfer  functions,  which 
are  shown  in  Figures  4-17  and  4-18  respectively.  Both  the  HH  and  VV  MTF  show  a 
null  about  the  range  direction,  and  both  have  gentle  humps  on  either  side  of  the  range 
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Figure  4-17.  Magnitude  Squared  of  HH  MTF  for  Runs  6-4,  7-2  and  7-4  Integrated  Over 
Wavenumber 
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direction.  The  hump  for  angles  less  than  90"  is  larger  than  the  hump  for  angles 
greater  than  90°  in  the  case  of  HH  polarization,  while  the  opposite  is  true  for  VV 
polarization. 

Comparing  Figures  4-16  and  4-18  with  Figures  4-15  and  4-16  respectively,  it  can 
be  seen  that  there  is  a  direct  correlation  between  the  angular  behavior  of  the  image 
spectra  and  the  angular  behavior  of  the  MTFs.  For  instance,  the  peak  of  the  HH 
image  spectra  for  Run  7-4  is  larger  than  that  for  Run  7-2  since  the  HH  MTF  is  larger 
for  azimuth  angles  less  than  90".  Similarily,  the  null  in  the  Run  6-4  VV  im.age 
spectrum  is  deeper  than  the  null  in  the  HH  image  spectrum  since  the  null  in  the  VV 
MTF  is  deeper  than  that  in  the  HH  MTF. 

The  asymmetrical  behavior  of  the  MlPs  as  a  function  of  azimuth  angle  is 
interesting  since  the  magnitude  of  the  individual  components  of  the  MTF  are 
symmetric  about  the  range  direction,  as  shown  in  Figure  4-19.  The  asymmetrical 
sum  of  the  component  MTFs  results  from  the  phasing  of  the  individual  components. 
Tests  showed  that  the  angular  behavior  of  the  MTF,  and  thus  the  angular  behavior  of 
the  image  spectrum,  could  be  significantly  altered  by  simply  varying  the  phasing  of 
the  individual  MTF  components. 

4.2.  WAVE  AMPLITUDE  ESTIMATION 

Having  gained  some  confidence  in  the  linear  imaging  model  as  well  as  the 
parameters  entering  this  model,  we  now  consider  the  inverse  problem  of  estimating 
wave  amplitudes  from  SAR  images.  The  general  procedure  we  followed  was  to 
Fourier  transform  the  image  intensity,  make  corrections  for  the  SAR  MTF  in  the 
wavenumber  domain  while  maintaining  the  phase  of  the  Fourier  transform,  and  then 
inverse  transform  to  obtain  an  estimate  of  the  surface  height  field.  However  some 
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Figure  4-19.  Magnitude  Squared  of  Component  MTFs  for  Run  6-4  Integrated  Over 
Wavenumber 
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preliminary  processing  steps  were  also  required,  as  described  in  the  following 
paragraph. 

The  received  SAR  intensity  (magnitude  squared)  image  I(x,y)  is  first  resampled 
from  the  slant  plane  to  the  ground  plane.  This  ground-sampled  data  is  then  passed 
through  a  median  filter  to  estimate  the  large  scale  trends  which  appear  in  the  imagery 
due  to  either  antenna  pattern  weighting  or  large  scale  scattering  changes  across  the 
image.  The  application  of  the  median  filter  is  simply  a  moving  box  across  the  image 
in  which  the  value  of  the  center  pixel  of  the  box  is  replaced  with  the  median  value  of 
the  surrounding  data.  For  this  analysis,  the  moving  box  had  an  extent  of  about  200 
meters  in  azimuth  and  200  meters  in  range.  The  images  used  for  this  analysis  cover 
approximately  3  kilometers  in  the  azimuth  direction  and  3  kilometers  in  the  range 
direction.  A  median  filter  was  used  as  opposed  to  a  boxcar  filter  because  the  median 
filter  retains  sharpness  of  features  (e.g.  the  centerline  wake)  whereas  the  boxcar  filter 
smooths  these  features  out.  The  unfiltered  image  is  then  divided  by  the  filtered  image 
on  a  pixel-by-pixel  basis  to  generate  a  "flattened”  image.  Clearly  this  process 
removes  low  frequency  information  from  the  image,  however,  the  range  of 
wavenumbers  occupied  in  the  Fourier  space  by  the  Kelvin  wake  is  far  above  those 
filtered.  This  filtered  image  is  then  normalized  to  have  unity  mean  and  the  FFT  is 
generated. 

Recall  from  Section  2.3  that  the  Kelvin  wake  occupies  a  narrow  hourglass  shape 
in  the  spectrum  whose  position  is  described  by  the  ship  speed  and  heading.  To 
estimate  the  Kelvin  wave  height,  the  spectral  energy  due  to  the  Kelvin  wake  needs  to 
be  extracted  so  that  the  inverse  transform  can  be  taken  to  generate  a  "waveheight 
image".  The  algorithm  for  extracting  the  Kelvin  spectral  energy  is  straightforward. 
The  energy  along  the  Kelvin  wake  is  summed  as  a  function  of  wavenumber  along  the 
wake  response.  Note  that  the  ship  velocity  displaces  the  hourglass  shape  closer  to  and 
further  away  from  the  DC  for  higher  and  lower  ship  speeds  respectively.  A  one-half 
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knot  width  surrounding  the  actual  ship  speed  will  sufficiently  contain  the  Kelvin 
energy  in  the  spectrum  (i.e.  for  a  ship  traveling  at  22  knots,  the  energy  is  summed 
over  a  wavenumber  range  corresponding  to  a  velocity  of  22  ±  .25  knots  at  the 
specified  heading).  This  energy  is  summed  separately  for  both  positive  and  negative 
propagation  angles,  (i.e.  both  arms  of  the  wake).  Once  the  energy  is  calculated  as  a 
function  of  wavenumber,  the  values  in  each  wavenumber  bin  are  normalized  for  the 
number  of  contributors  to  that  bin  so  that  an  average  energy  level  is  calculated  for  the 
spectrum.  In  a  likewise  manner,  the  Kelvin  waveheight  spectrum  can  be  calculated  as 
a  function  of  spectral  angle  (/»=tan  ‘(Kj,/KJ  as  well. 

A  correction  was  also  made  to  remove  any  spectral  bias  present  that  might 
contribute  energy  to  the  Kelvin  spectral  response.  If  the  spectra  of  interest  were  from 
wavefields  with  a  distributed  range  of  wavenumber  components  one  t.'^ipht  refer  to  the 
algorithms  described  earlier  to  estimate  the  spectral  bias.  Because  the  Kelvin  wake 
has  a  very  narrow  spectral  response,  the  approach  for  estimating  the  Kelvin 
waveheight  spectrum  can  be  extended  to  estimate  the  ambient  spectrum.  The 
approach  taken  is  to  generate  an  ambient  spectrum  over  the  Kelvin  hourglass  shape, 
but  at  speeds  which  are  the  known  ship  speed  ±  2  knots.  It  is  assumed  that  the 
ambient  spectrum  will  not  change  significantly  over  the  4  knot  width  used  to  estimate 
the  different  waveheight  spectra.  This  assumption  has  been  validated  for  the  data 
examined.  Note  that  for  different  ship  speeds,  the  range  of  wavenumbers  covered  by 
the  hourglass  shape  vary.  To  allow  for  a  direct  comparison  between  the  different 
waveheight  estimates  (both  the  wake  and  the  background  estimates)  it  is  necessary  to 
normalize  the  different  spectra  to  have  equivalent  wavenumbers.  This  is 
accomplished  by  scaling  the  wavenumbers  for  each  spectrum  by 

k’  =  kV.^/g.  (19) 
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This  normalization  allows  direct  comparison  to  be  made  between  the  different 
waveheight  spectra.  If  the  data  is  generated  as  a  function  of  spectral  angle,  no 
normalization  is  necessary.  Once  the  background  waveheight  spectra  have  been 
generated,  they  are  averaged  together  to  generate  a  single  spectral  background 
estimate.  This  background  estimate  is  smoothed  and  subtracted  from  the  slightly 
smoothed  Kelvin  wake  waveheight  spectrum  to  generate  a  final  waveheight  spectrum 
estimate.  An  example  of  the  waveheight  spectra  as  a  function  of  angle  is  shown  in 
Figure  4-20  where  the  energy  calculated  over  the  Kelvin  response  is  shown  along  with 
the  ambient  spectrum  estimate.  The  same  spectrum  plotted  as  a  function  of 
normalized  wavenumber,  k’,  is  shown  in  Figure  4-21. 

To  generate  the  actual  2-dimensional  waveheight  image  it  is  necessary  to  use  the 
two  dimensional  complex  FFT.  The  2-D  FFT  is  bandpassed  over  the  Kelvin 
hourglass  shape  with  a  spectral  width  corresponding  to  the  ±  0.25  knots  described 
earlier.  The  spectral  background  estimate  is  then  subtracted  over  this  band  to 
generate  the  final  two-dimensional  FFT.  Because  the  FFT  is  complex  and  contains 
essential  phase  information,  the  subtraction  is  performed  using  the  intensity  of  the 
data.  This  data  is  then  converted  back  to  complex  form  by  taking  the  square  root  of 
the  intensity  and  applying  the  phase  of  the  original  FFT  to  each  pixel. 

At  this  point  the  correction  for  the  SAR  MTF  is  finally  made.  Because  the  MTF 
is  not  symmetric  it  is  important  to  know  which  quadrants  of  the  modeled  transfer 
function  are  effecting  the  actual  spectrum.  Each  quadrant  of  the  MTF  represents  a 
particular  wave  direction.  A  visual  examination  of  the  data  indicates  which  quadiants 
of  the  MTF  should  be  used  (based  on  the  ship  heading)  and  a  Hermitian  version  of 
the  MTF  is  generated  by  mirroring  the  chosen  quadrants  and  conjugating.  This  MTF 
is  coherently  divided  out  of  the  existing  wake  spectrum  to  generate  a  final  estimate  of 
the  2-D  Kelvin  waveheight  spectrum.  The  data  is  inverse  transformed  after  zeroing 
out  the  DC  term  and  the  real  part  of  the  output  represents  the  actual  waveheight. 
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Figure  4-20.  Image  Spectral  Densities  Along  Kelvin  Wake,  Plotted  Versus  Wake  Angle 
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Scans  are  generated  over  the  walce  feature  to  examine  the  waveheight  (which  is  in 
units  of  meters)  over  the  cusps  of  both  wake  arms. 

4.2.1.  Ambient  Wave  Results 

A  measurement  of  the  ambient  surface  waveheight  was  made  using  the  described 
algorithm  to  verify  its  capability.  The  ambient  data  for  Run  7-4  was  chosen  as  a  test 
case  because  a  relatively  large  wave  field  could  be  extracted  without  interference  from 
the  Kelvin  wake.  The  data  set  for  this  test  was  extracted  at  a  61  degree  incidence 
angle  with  an  R/V  ratio  of  42  seconds.  The  spatial  extent  of  this  data  is  1.5  km  in 
azimuth  and  2.5  km  in  ground  range. 

The  normalized  image  intensity  was  generated  for  both  X-band  and  L-band  data 
as  described  above.  The  Fourier  transform  was  generated  and  a  mean  background 
level  was  removed  from  the  spectra  to  compensate  for  the  speckle  contribution.  This 
background  value  was  calculated  as  the  mean  spectral  value  in  a  homogeneous  area 
near  the  peak  spectral  response.  Naturally,  the  speckle  level  will  be  over-estimated  at 
the  higher  wavenumbers  and  under-estimated  at  the  lower  wavenumbers  due  to  the 
system  impulse  response.  This  will  not  effect  the  measurement  significantly, 
however,  because  the  bandpass  region  chosen  is  taken  over  the  spectral  peak  and  has 
a  small  wavenumber  extent.  The  estimated  SAR  MTF  effects  are  removed  using  the 
previously  described  model.  The  data  is  symmetrized  and  a  bandpass  region  is 
chosen  which  surrounds  the  peak  ambient  spectral  response.  The  peak  wavelength 
was  measured  to  be  88  m  at  an  angle  of  65  degrees  from  azimuth.  The  bandpass  for 
this  data  was  taken  to  be  a  small  spectral  box  surrounding  the  peak  spectral  response 
with  an  azimuth  wavenumber  extent  of  0.057  rads/m  and  a  range  wavenumber  extent 
of  0.038  rads/m.  The  data  was  then  inverse  transformed  to  form  the  waveheight 
image.  The  original  image  intensities  are  shown  in  figures  4-22  and  4-23  with  the 


65 


■^RIM 


Figure  4-22.  (a)  Image  Subset,  (b)  Image  Spectrum,  and  (c)  Wave  Height  Image  for  Run 
7-4.  X-VV 
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resulting  :  ctra  and  waveheight  images.  The  RMS  value  of  the  waveheight  is 
calculated  over  the  entire  image  to  be  0.222  m  at  X-band  and  0.200  m  at  L-band. 
These  vdues  compare  well  with  the  value  of  0.23  m  measured  in  situ  (Gasparovic, 
1989).  The  waveheight  estimated  from  the  X-band  image  data  compares  well  to  the 
waveheight  measured  in  situ  although  the  L-band  estimate  is  slightly  lower  than  the 
in  situ  values.  Scans  through  this  wavefield  are  shown  in  Figures  4-24  and  4-25.  A 
possible  cause  for  the  reduced  waveheight  may  be  a  misfocussing  of  the  waves  due  to 
larger  integration  times  at  L-band. 

An  important  parameter  in  the  MTF  model  is  the  standard  deviation  of  the  r.m.s. 
surface  velocity,  which  is  a  function  of  windspeed  and  affects  the  azimuthal  cutoff 
of  the  spectrum.  This  value  is  unknown  for  this  data  set  and  as  a  result,  waveheight 
measurements  were  generated  for  a  range  of  values  at  X-band  to  determine  how 
sensitive  the  measured  waveheight  is  to  o,,.  Figure  4-26  shows  how  the  MTF 
amplitude  changes  as  a  function  of  a,,.  Note  that  the  azimuth  bandwidth  (azimuth  is 
horizontal)  decreases  quickly  with  increasing  o;  which  indicates  that  as  the  wind  speed 
increases,  the  shorter-scale  azimuth-travelling  waves  cannot  be  resolved.  Figure  4-27 
contains  the  measured  waveheight  as  a  function  of  o;.  This  curve  indicates  that  the 
measured  long-scale  waveheight  is  somewhat  insensitive  to  r.m.s.  velocity  variations. 
This  result  is  obtained  because  the  waves  are  close  to  range  traveling  in  direction.  As 
the  wavefronts  become  more  azimuth  traveling,  the  azimuth  cut-off  would 
significantly  effect  the  ability  to  resolve  the  shorter  waves. 

4.2.2.  Kelvin  Wake  Results 

The  algorithm  was  applied  to  three  data  sets  which  contain  Kelvin  wakes.  These 
data  sets  all  contain  destroyer  (DDG)  wakes  with  two  containing  azimuth  traveling 
ships  (one  being  near-shore,  one  off-shore)  and  the  last,  a  range-traveling  case. 
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Figure  4-25.  Ambient  Waveheight  Scan  for  Run  7-4,  L-VV 
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Figure  4*27.  Estimated  Waveheight  Versus  Assumed  Value  of  r.m.s.  Surface  Velocity 
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Imagery  of  these  areas  are  shown  in  Figures  3-3  through  3-6.  Subsets  of  this  imagery 
taken  around  the  wake  feature  were  used  for  analysis.  The  algorithm  described  above 
was  applied  to  these  data  sets  and  the  results  follow. 

Figure  4-28  contains  the  intensity  image,  resulting  spectrum,  and  calculated 
waveheight  image  for  DDG  Run  4-1  (X-Band  VV)  which  is  the  near-shore  wake 
image.  The  hourglass  shape  of  the  Kelvin  waxe  spectral  response  is  clearly  visisble. 
Because  the  value  of  is  unknown,  waveheight  estimates  are  calculated  for  two 
values  of  which  bracket  a  reasonable  value  for  this  parameter,  ff„=0.2  m/s  and 
<Tv=0.4  m/s.  Figure  4-29  contains  waveheight  measurements  for  DDG  Run  4-1  (X- 
VV)  for  both  the  Port  and  Starboard  Arms  along  the  cusp  of  the  wake  at  both 
values.  There  is  a  significant  difference  in  the  measured  waveheight  for  these 
different  values  of  (t„.  At  0^  =  .!  m/s  the  peak  waveheight  for  both  arms  of  the  wake 
is  about  0.2  meters  while  this  value  is  about  .5  meters  for  (t,=.4  m/s.  Note  the 
amplitude  falls  off  with  increased  distance  behind  the  wake  along  this  cut.  The 
corresponding  L-Band  imagery  and  scans  are  shown  in  figures  4-30  and  4-31 
respectively.  The  measured  L-band  waveheight  is  smaller  than  that  of  X-band  with 
the  peak  waveheights  beinj  about  0.1  meters  for  a^-.2  m/s  and  about  0.4  meters  for 
ffv=.4  m/s.  Note  that  for  the  L-band  data,  the  values  are  not  the  same  for  both  the 
port  and  st'  .'  oard  arms.  This  asymmetr}-  may  be  due  to  a  difference  in  surface 
conditions  (sea  state)  between  the  two  areas. 

Figure  4-32  contains  the  intensity  image,  resulting  spectrum,  and  calculated 
waveheight  image  for  DDG  Run  7-4  (X-Band  VV)  whi^h  is  the  open-water  azimuth 
traveling  wake  image.  The  Kelvin  wake  is  less  pronounced  than  in  Run  4-1  but  is 
still  clearly  visible.  Scans  of  the  X-band  waveheight  are  shown  in  figure  4-33  for  the 
two  values  of  a^.  The  measurements  show  a  peak  waveheight  of  about  0.15  meters 
for  ay=.2  m/s  and  0.2  meters  for  (t„  =  .4  m/s.  The  difference  between  these  values 
and  those  for  Run  4-1  may  be  due  to  sea  state  effects.  The  corresponding  L-band 
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Figure  4-29.  Kelvin  Waveheight  Estimates  for  Cusp  Waves  Along  Port  and 
Starboard  Arms  for  Pass  4-1  Using  ov  =  0.2  m/s  and  0.4  m/s 
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Figure  4-30.  (a)  Image  Subset,  (b)  Image  Spectrum,  and  (c)  Kelvin  Waveheight  Image  for 
Run4-1,L-VV 
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Figure  4-31 .  Kelvin  Waveheight  Estimates  for  Cusp  Wave;;  Along  Port  and 
Starboard  Arms  for  Pass  4-1  Using  ov  =  0.2  m's  and  0.4  m/s 
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Figure  4-32.  (a)  Image  Subset,  (b)  Image  Spectrum,  and  (c)  Kelvin  Waveheight  Image  for 
Run  7-4,  X-VV 
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Figure  4*33.  Kelvin  Waveheight  Estimates  for  Cusp  Waves  Along  Port  and 
Starboard  Arms  for  Pass  7-4  Using  ov  =  0.2  m/s  and  0.4  m/s 
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results  are  shown  in  figures  4-34  and  4-35  which  show  the  same  characteristics  as  the 
X-band  data. 

Figures  4-36  and  4-38  contain  the  intensity  image,  resulting  spectrum,  and 
calculated  waveheight  image  for  DDG  Run  7-2  (X-Band  VV  and  L-Band  VV 
respectively)  which  is  the  open-water  range-traveling  wake  image.  Scans  of  the 
calculated  waveheight  are  shown  in  figures  4-37  and  4-39  for  both  the  X-band  and 
L-band  data  sets  for  <Tv=.2  m/s  only.  These  figures  show  the  peak  waveheight  to  be 
about  .15  m  for  X-band  data  and  about  .10  m  for  L-band  data  along  the  cusp  although 
waveheights  as  large  as  0.3  m  are  measured  for  the  stem  wave  component  directly 
behind  the  ship.  Because  of  the  SAR  MTF  effects,  the  predominant  Kelvin  wake 
component  is  the  transverse  component  which  leads  to  relatively  more  wake  energy  in 
the  stem  wave  component. 
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Figure  4-34.  (a)  Image  Subset,  (b)  Image  Spectrum,  and  (c)  Kelvin  Waveheight  Image  for 
Run  7-4,  L-VV 
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Figure  4-35.  Kelvin  Waveheight  Estimates  for  Cusp  Waves  Along  Port  and 
Starboard  Arms  for  Pass  7-4  Using  ov  =  0.2  m/s  and  0.4  m/s 
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Figure  4-36.  (a)  Image  Subset,  (b)  Image  Spectrum,  and  (c)  Kelvin  Waveheight  image  for 
Run  7-2,  X-VV 
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Figure  4-37.  Kelvin  Waveheight  Estimates  for  Cusp  Waves  Along  Port  and 

Starboard  Arms  for  Run  7-2,  X-W  Using  ov  =  0.2  m/s  and  0.4  m/s 
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Figure  4-38.  (a)  Image  Subset,  (b)  Image  Spectrum,  and  (c)  Kelvin  Waveheight  Image  for 
Run7-2,  L-VV 
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Figure  4-39.  Kelvin  Waveheight  Estimates  for  Cusp  Waves  Along  Port  and 
Starboard  Arms  for  Pass  7-2  Using  ov  =  0.2  m/s  and  0.4  m/s 
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5.0  SUMMARY  AND  CONCLUSIONS 

The  dependence  of  the  SAR  modulation  transfer  function  on  a/imuth  angle  and 
polarization  was  investigated  by  examining  ocean  wave  image  spectra  extracted  from 
collected  imagery  and  by  modeling  the  corresponding  image  spectra.  Good  agreement 
was  observed  between  the  empirical  and  modeled  spectra  although  some  differences 
were  observed  near  the  range  wavenumber  axis,  probably  due  to  errors  in  the 
hydrodynamic  modulation  transfer  function.  The  close  correspondence  suggests  that 
the  model  of  the  S.\R  modulation  transfer  function  summarized  in  section  2.1  is 
reasonably  accurate  for  the  conditions  under  which  this  data  was  collected. 

The  success  of  the  modeling  effort  also  encouraged  us  tr  attempt  an  estimation  of 
the  wave  height  field  for  both  ambient  and  Kelvin  wake  waves.  The  ambient  wave 
height  estimates  were  in  good  agreement  with  measured  wave  heights.  Altnough  no 
measurements  of  the  DDG  Kelvin  wave  heights  were  available  for  compai'son,  the 
values  estimated  from  the  SAR  imagery  were  felt  to  be  reasonable.  However,  it 
should  be  noted  that  the  estimated  values  were  in  some  cases  quite  sensitive  to  the 
assumed  value  of  the  surfo„c  radial  velocity  variance.  Since  this  parameter  is  not 
well  known,  this  leads  to  a  significant  uncertainty  in  the  wave  height  estimates  for 
some  viewing  geometries. 
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